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Abstract. By manipulating the spherical Jeans equation, Wolf et al. (2010) show that the mass 
enclosed within the 3D deprojected half-light radius r 1/2 can be determined with only mild 
assumptions about the spatial variation of the stellar velocity dispersion anisotropy as long as 
the projected velocity dispersion profile is fairly flat near the half-light radius, as is typically 
observed. They find M 1/2 = 3G" 1 (afj) r 1/2 ~ 4G" 1 (of oa }R c , where (cr^J is the luminosity- 
weighted square of the line-of-sight velocity dispersion and R c is the 2D projected half-light 
f**> [■ radius. This finding can be used to show that all of the Milky Way dwarf spheroidal galaxies 
(MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10 9 Mq , 
assuming a ACDM cosmology. In addition, the dynamical I-band mass-to-light ratio T vs. 
„H ■ M 1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near 
Y; ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to Y' ~ 3,200 
' for ultra-faint dSphs, and a more shallow rise to Y\ /2 — 800 for galaxy cluster spheroids. 

' Keywords, galaxies: formation, kinematics and dynamics, Local Group, 
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1. Introduction 

Mass determinations for dispersion-supported galaxies based on only line-of-sight ve- 



locity measurements suffer from an uncertainty associated with not knowing the intrinsic 
3D velocity dispersion. The difference between tangential and radial velocity disper- 
sions is quantified by the stellar velocity dispersion anisotropy, /?. Many questions in 
galaxy formation are affected by our ignorance of /3, including the ability to quantify the 



■ dark matter content in the outer parts of elliptical galaxies Romanowsky et al. (2003) 
| Dekel et al. (2005)] , to measure the mass profile of the Milky Way from stellar halo kine 



£N) . matics Battaglia et al. (2005)| |Dehnen et al. (2006)] , and to infer accurate mass distri 



butions in dwarf spheroidal galaxies (dSphs) |Gilmore et al. (2007) Strigari et al. (2007)] . 



. Wolf et al. (2010) (hereafter W10) used the spherical Jeans equation to show that for 
each dispersion-supported galaxy, there exists one radius within which the integrated 
mass as inferred from the line-of-sight velocity dispersion is largely insensitive to (3, 
and that for a wide range of stellar distributions, this radius is approximately the 3D 
deprojected half-light radius r 1/2 : 

M^EMfr.J^G-'^Jv, (1.1) 



5_j 

~4G- 1 (a i 2 oa )R 



930 (Ml,) 1 M 



\ km 2 s 



pc 



(of ) is the luminosity- weighted square of the line-of-sight velocity dispersion. In the 
second line we have used R o ~ (3/4)r 1/2 for the 2D projected half-light radius. This 
approximation is accurate to better than 2% for exponential, Gaussian, King, Plummer, 
and Sersic profiles (see Appendix B of W10 for useful fitting formulae). 
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2. The Spherical Jeans Equation 

Given the relative weakness of the scalar virial theorem as a mass estimator (see Section 
2 of W10), we turn to the spherical Jeans equation. It relates the tracer velocity dispersion 
and tracer number density n*(r) of a spherically symmetric, dispersion-supported, colli- 
sionless stationary system to its total gravitating potential $(r), under the assumption 
of dynamical equilibrium with no streaming motions: 

^dg _ d(n«<r?) | ^n,j r 2 
dr dr r 

<r r (r) is the radial velocity dispersion of the stars/tracers and (3{r) = 1 — of/cr 2 measures 
the velocity anisotropy, where the tangential velocity dispersion <j t — erg = a^. It is 
informative to rewrite the implied total mass profile as 

M(r) = T -^ ( 7 * + 7a -2/3), (2.2) 

where 7* = — dlnn*/dlnr and 7 CT = — dlncr 2 /dlnr. Given line-of-sight kinematics, the 
only term on the right-hand side of Equation l2.2l that can be determined by observations is 
7*, which follows from the projected surface brightness profile under an assumption about 
how it is related to the projected stellar number density £*(ii). Via an Abel inversion 
n* is mapped in a one-to-one manner with the spherically deprojected observed surface 
brightness profile by assuming that n± traces the light density. 

Line-of-sight kinematic data provides the projected velocity dispersion profile a los (R). 



As first shown by Binney & Mamon (1982) in order to use the Jeans equation one must 
relate a, , to a r : 



/•oo 

E*(7. a (R) = / n*a 2 r {r) 
Jr 2 
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1 - 5-0(r) 



dr 2 



Vr 2 - R 2 



(2.3) 



It is clear then that a significant degeneracy associated with using the observed £*(i?) 
and cr loa (i?) profiles exists in trying to determine an underlying mass profile M(r) at any 
radius, as uncertainties in j3 will affect both the relationship between M(r) and oy in 
Equation 12.21 and the mapping between a r and cr los in Equation 12.31 

The technique of W10 for handling the j3 degeneracy in order to provide a fair rep- 
resentation of the allowed mass profile given a set of observables is to consider general 
parameterizations for M(r) and j3{r) and then to undertake a maximum likelihood anal- 
ysis to constrain all possible parameter combinations. By using such a strategy, W10 
derive meaningful mass likelihoods for a number of dispersion-supported galaxies with 
line-of-sight velocity data sets. 

The stellar velocity dispersion anisotropy can be modeled as a three-parameter function 

/9(r) = (A»- #0-3^7-5- + A>, (2-4) 

and the total mass density distribution can be described using the six-parameter function 

p s e" r/r cut 

Ptot(r) = (r/r.)T[l + (r/r a ) a ]('-T)/«' (2 ' 5) 

For their marginalization, W10 adopt uniform priors over the following ranges: log 10 (r 1/2 /5) < 
log 10 Og < log 10 (r u J; -10 < /3oo < 0.91; -10 < ft, < 0.91; log 10 (r 1/2 /5) < log 10 (r s ) < 
log 10 (2 r ); 0<7<2;3<^<5; and 0.5 < a < 3, where r lim is the truncation radius 
for the stellar density. The variable r cut allows the dark matter halo profile to truncate 
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Figure 1. Figure 1 from W10. See the text for details. Right: The green dot-dashed line repre- 
sents the contribution of mass from the stars, assuming a stellar V-band mass-to-light ratio of 
3 M /L . 



beyond the stellar extent and the uniform prior log 10 (r lim ) < log 10 (r cut ) < log 10 (r h h ) 
is adopted. For distant galaxies W10 use r hi h = 10 r lim and for satellite galaxies of the 
Milky Way they set r high equal to the Roche limit for a 10 9 M Q point mass. In practice, 
this allowance for r cut is not important when focusing on integrated masses within the 
stellar radius. 

W10 apply their marginalization procedure to resolved kinematic data for MW dSphs, 
MW globular clusters, and elliptical galaxies. Since MW dSphs and globular clusters are 
close enough for individual stars to be resolved, the joint probability of obtaining each 
observed stellar velocity given its observational error and the predicted line-of-sight ve- 
locity dispersion from Equations 12.11 and 12.31 is considered. In modeling the line-of-sight 
velocity distribution for any system, the observed distribution is a convolution of the 
intrinsic velocity distribution, arising from the distribution function, and the measure- 
ment uncertainty from each individual star. Given that line-of-sight velocity distributions 
of dispersion-supported systems are often well-described by a Gaussian, the probability 
of obtaining a set of line-of-sight velocities Y given a set of model parameters ^ is 
described by the likelihood 



P{r\~#) = n 



: exp 



1 {% ~ vf 



2 <s 



(2.6) 



The product is over the set of N stars, where v is the average velocity of the galaxy. As 
expected, the total error at a projected position is a sum in quadrature of the theoretical 
intrinsic dispersion, Cth,i(^), and the measurement error e.;. The posterior probability 
distribution for the mass at any radius can be generated by multiplying the likelihood 
by the prior distribution for each of the nine /9tot(?~) and /3(r) parameters as well as 
the observationally derived parameters and associated errors that yield n+(r) for each 
galaxy, including uncertainties in distance. We then integrate over all model parameters, 
including v, to derive a likelihood for mass. Following |Martinez et al. (2009) a Markov 



Chain Monte Carlo technique is used in order to perform the required ten to twelve 
dimensional integral. 
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For elliptical galaxies that are located too far for individual stellar spectra to be ob- 
tained, the resolved dispersion profiles are analyzed with the likelihood 



f 



1 (g» - Qth^f 

2 e? 



where the product is over the set of N dispersion measurements 
error of each measurement. 



(2-7) 

and ei is the reported 



3. Minimizing the Anisotropy Degeneracy 

As discussed in W10, the degeneracy between the anisotropy and the integrated mass 
will be minimized at an intermediate radius within the stellar distribution. Such an 
expectation follows from considering the relationship between ay and cr los . 

At the observed center of a spherical, dispersion-supported galaxy (R = 0), line-of- 
sight velocities project onto the radial component with <r los ~ ay, while at the edge of 
the system (R = r Um ), line-of-sight velocities project onto the tangential component with 
<r los ~ at- As an example, consider an intrinsically isotropic galaxy (/3 = 0). If this system 
is analyzed using line-of-sight kinematics under the false assumption that ay < at (ft < 0) 
at all radii, then the total velocity dispersion at r ~ would be overestimated while the 
total velocity dispersion at r ~ r lim would be underestimated. Conversely, if one were to 
analyze the same galaxy under the assumption that ay > a t (ft > 0) at all radii, then the 
total velocity dispersion would be underestimated near the center and overestimated near 
the galaxy edge. Thus, some intermediate radius should exist where attempting to infer 
the enclosed mass from only line-of-sight velocities is minimally affected by the unknown 
value of ft. 

These qualitative expectations are quantitatively displayed in Figure [1] where inferred 
mass profiles for the Carina dSph galaxy for several choices of constant ft are shown. The 
right-hand panel shows the same data analyzed using the complete likelihood analysis, 
where the fairly general ft(r) profile presented in Equation 12.41 was marginalized over. 
The dataset is discussed in W10, where the average velocity error is ~ 3 km s _1 . Each line 
in the left panel of Figure Q] shows the median likelihood of the cumulative mass value at 
each radius for the value of ft indicated. The 3D half-light radius and the limiting stellar 
radius are marked for reference. As expected, forcing ft < produces a systematically 
higher (lower) mass at a small (large) radius compared to ft > 0. Thus, this requires that 
every pair of M(r) profiles analyzed with different assumptions about ft cross at some 
intermediate radius. Somewhat remarkable is that there exists a radius where every pair 
approximately intersects. We see that this radius is very close to the deprojected 3D 
half-light radius r 1/2 . The right-hand panel in Figure [T] shows the full mass likelihood as 
a function of radius, with the shaded bands illustrating the 68% and 95% likelihood con- 
tours. The likelihood contour also pinches near r 1/2 , as the data preferentially constrain 
this mass value. 

By examining each of the well-sampled dSph kinematic data sets in more detail, W10 
finds that the measurement errors, rather than the ft uncertainty, always dominate the 
errors on the mass near r 1/2 , while the mass errors at both smaller and larger radii are 
dominated by the ft uncertainty (and thus are less affected by measurement error). To 
analytically describe this effect, let us briefly examine the Jeans equation in the context 
of observables. 

Consider a dispersion-supported stellar system where S* (R) and ct 1ob (R) are deter- 
mined accurately by observations, such that any viable solution will keep the quantity 
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er 2 s (i?) fixed to within allowable errors. With this in mind, W10 show that Equa- 
tion [273] can be rewritten in a form that is invertible, isolating into a kernel the integral's 
i?-dependence: 



ii 1 



n*<7 2 (r) 



R 2 



dr 2 



V V 2 



n*ol (1 - /3)r 2 + /3(r 2 - R 2 ) 



i? 2 



i? 2 



dr 2 



(3.1) 
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fl 2 
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dr 2 
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2 V^ 2 



R 2 



2r 



df 2 



dr 2 



VT 



i? 2 



where an integration by parts was utilized to achieve the third equality. Note that the 
second term on the third line must be zero under the physically-motivated assumption 
that the combination Pn+a 2 falls faster than r _1 at large r.0 

Because Equation 13.11 is invertible, the fact that the left-hand side is an 
observed quantity and independent of /3 implies that the term in brackets 
must be well determined regardless of a chosen f}. This allows one to equate the 
isotropic integrand with an arbitrary anisotropic integrand: 



r l,a=o 



n*a r 2 [l-/3(r)] 



/3n*cr 2 df 



(3.2) 



After taking a derivative with respect to lnr and subtracting Eq uation 12 . 1 1 the following 
result is obtained 



M(r; (3) - M(r; 0) = (Mj^M ^ + ^ + ^ _ 3) 



(3.3) 



We remind the reader that 7* = — dlnn*/dlnr and 7^ = — dln<r 2 /dlnr. Likewise, 
7/3 = — dln./3/dlnr = —f3'//3, where ' denotes a derivative with respect to lnr. The 
possibility of a radius r oq is revealed by Equation l3.31 where the term in parentheses goes 
to zero, such that the enclosed mass M(r oq ) is minimally affected by not knowing /3(r): 



7*(r cq ) = 3 - 7<x(r cq ) ~ J/3 (r oq ) 



(3.4) 



W10 go through additional detailed arguments to justify that r cq ~ r 1/2 for systems 
with relatively flat observed dispersion profiles. We will utilize this finding (Equation ll.ll) 
to perform tests of galaxy formation. 



4. Milky Way dwarf spheroidal galaxies 

As an example of the utility of M determinations, Figure [5] presents M vs. r 1/2 for 
MW dSph galaxies. Relevant parameters for each of the galaxies are provided in Table 1 
of W10. The symbol types labeled on the plot correspond to three luminosity bins that 
span almost five orders of magnitude in luminosity. It is interesting to compare the data 

f Appendix A of W10 shows how an Abel inversion can be used to solve for oy(r) and M(r) 
in terms of directly observable quantities. 
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Figure 2. Figure 3 from W10. The half-light masses of the Milky Way dSphs plotted against 

r 1/2 . See the text for descriptions. 



points in Figure [2] to the integrated mass profile M(r) predicted for ACDM dark matter 
field halos of a given Mhaio mass, which is defined as the halo mass corresponding to an 
overdensity of 200 compared to the critical density. In the limit that dark matter halo 
mass profiles M(r) map in a one-to-one way with their Mhaio mass, then the points on 
this figure may be used to estimate an associated halo mass for each galaxy. 

The solid line in the left panel of Figure[2]shows the mass profile for a NFW dark matter 
halo |Navarro et al. (1997)] at z = with a halo mass Mhaio = 3 x 10 9 M©. The median 
concentration (c = 11) predicted by the mass-concentration model |Bullock et al. (2001)| 
updated for WMAP5 ACDM parameters [Maccio et al. (2008)] is used. The dashed lines 
indicate the expected 68% scatter about the median concentration at this mass. The 
dot-dashed line shows the expected M(r) profile for the same mass halo at z — 3 (corre- 
sponding to a concentration of c — 4) , which provides an idea of the spread that would 
result from the scatter in infall times. An interesting result is that each MW dSph is con- 



sistent with inhabiting a dark matter halo of mass ~3x 10 M© Strigari et al. (2008) 

The right panel in Figure [5] shows the same data plotted with the median mass profiles 
for different halo masses. Clearly, the MW dSphs are also consistent with populating 
dark matter halos of a wide range in Mhaio above ~ 3 x 10 8 M©. This result informs of 
a very stringent limit on the fraction of the baryons converted to stars in these systems. 
More importantly, there is no systematic correlation between luminosity and the Mhaio 
mass profile that each dSph most closely intersect^]. The ultra-faint dSph population 
(circles) with Ly < 10 4 L© is equally likely to be associated with the more massive 
dark matter halos as are classical dSphs that are more than three orders of magnitude 
brighter (squares). A simple interpretation of the right-hand panel of Figure [5] shows that 
the two least luminous satellites (which also have the smallest M. .. and r, .„ values) are 

\ 1/ 2 1/2 i—i ' 

associated with halos that are more massive than any of the classical MW dSphsyj. This 
behavior is difficult to explain in models constructed to reproduce the Milky Way satellite 
population, which typically predict a trend between dSph luminosity and halo infall mass. 



t There are hints by Kalirai et al. (2010) 



X Martinez et al. (20lTTJ[ |Minor et al. (20 



that thi s result does not hol d for the M31 dSphs. 



0) and Simon et al. (2010) have explored the ef- 
fects of the inflation of the observed velocity dispersion due to binaries, but find that the faintest 
system, Segue 1, is still extremely dark-matter dominated. 
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It is possible that a new scale in galaxy formation exists Strigari et al. (2008)| or that 
there is a systematic bias that makes it difficult to detect low luminosity galaxies which 
sit within low mass halos |Bullock et al. (2010)] . 



5. Global population of dispersion-supported systems 

Figure [3] examines the relationship between the half-light mass M 1/2 and the half-light 
luminosity L 1/2 = 0.5 Li for the full range of dispersion-supported stellar systems in the 
universe: globular clusters, dSphs, dwarf ellipticals, ellipticals, brightest cluster galaxies, 
and extended cluster spheroids. 

There are several noteworthy aspects to Figure [3l each highlighted in a different way 
in the three panels. In the middle and right panels the half-light mass-to-light ratios of 
spheroidal galaxies in the universe demonstrate a minimum at T* /2 ~ 2 — 4 that spans 
a broad range of luminosities Lj ~ 10 85 ~ la5 L Q and masses M 1/2 ~ 10 9-11 M Q . It is 
interesting to note the offset in the average mass-to-light ratios between ellipticals and 
globular clusters, which may suggest that even within r 1/2 , dark matter may constitute 
the majority of the mass content of L* elliptical galaxies. Nevertheless, it seems that 
dark matter plays a dominant dynamical role (T* > 5) within r 1/2 in only the most 
extreme systems. The dramatic increase in half-light mass-to-light ratios at both larger 
and smaller luminosity and mass scales is indicative of a decrease in galaxy formation 
efficiency in the smallest^] and largest dark matter halos. It is worth mentioning that 
if ACDM is to explain the luminosity function of galaxies a similar trend in the Mhaio 
vs. L relationship must exist. While a different mass variable is presented in Figure [3l 
the resemblance between the two relationships is striking, and generally encouraging for 
ACDM. 

One may gain some qualitative insight into the physical processes that drive galaxy 
formation inefficiency in bright vs. faint systems by considering the L 1/2 vs. M 1/2 relation 
(left panel) in more detail. There exist three distinct power-law regimes M 1/2 cx L^ /2 
with H > 1, H ~ I, and K < 1 as mass decreases. Over the broad middle range of galaxy 
masses, M 1/2 ~ 10 9 " 11 M Q , mass and light track each other closely with K ~ 1, while 
faint galaxies obey K ~ 1/2, and bright elliptical galaxies are better described by H ~ 4/3 
transitioning to N > 1 for the most luminous cluster spheroids. One may interpret the 
transition from H > I in bright galaxies to N < I in faint galaxies as a transition 
between luminosity-suppressed galaxy formation to mass-suppressed galaxy formation. 
That is, for faint galaxies (H < I), there does not seem to be a low-luminosity threshold 
in galaxy formation, but rather behavior closer to a threshold (minimum) mass with 
variable luminosity. For brighter spheroids with K > 1, the increased mass-to- light ratios 
are driven more by increasing the mass at fixed luminosity, suggestive of a maximum 
luminosity scale (W10). 

Going a step further, Figure [3] provides a useful empirical benchmark against which 
theoretical models can compare. It will be challenging to explain how two of the least 
luminous MW dSphs have the highest mass-to-light ratios T 1 — 3, 200 of any the 

1/2 

collapsed structures shown, including intra-cluster light spheroids, which reach values of 
Y* /2 ~ 800. Not only are the ultra-faint dSphs the most dark matter dominated objects 



f These findings were expanded upon by |Tollerud et al. (2010) to show that all systems which 
sit embedded within dark matter halos lie along a one dimensional relation within three dimen- 
sion al space. 

£ |McGaugh fc Wolf (2010)] recently showed that correlations exist when one explores addi- 
tional observable properties of dSphs. 
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Figure 3. See the text for descriptions. The symbols are linked to galaxy types as follows: 
Milky Way dSphs (squares, diamonds, circles), galactic globular clusters (stars), dwarf ellipticals 
(triangles), ellipticals (inverted triangles), brightest cluster galaxies (plus signs), and cluster 
spheroids (asterisks). See Figure 4 of W10 for references. 

known, given that have even lower baryon-to-dark matter fractions /& ~ f^/Q^m J$ 10 -3 
than galaxy clusters /& ~ 0.1, W10 points out that ultra- faint dSphs also have higher 
mass-to- visible light ratios within their stellar extents than even the well-studied galaxy 
cluster spheroids. 
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